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Abstract. We consider the self-adjoint operator H — Ho + V , where Ho is the free semi- 
classical Dirac operator on R^. We suppose that the smooth matrix- valued potential V — 
0{{x)~^), S > 0, has an analytic continuation in a complex sector outside a compact. We 
define the resonances as the eigenvalues of the non-selfadjoint operator obtained from the 
Dirac operator H hy a complex distortions of K^. We establish an upper bound 0{h~^) 
for the number of resonances in any compact domain. For 5 > 3, a representation of the 
derivative of the spectral shift function ^(A, h) related to the semi-classical resonances of H 
and a local trace formula are obtained. In particular, if V is an electro-magnetic potential, 
we deduce a Weyl-type asymptotic of the spectral shift function. As a by-product, we obtain 
an upper bound 0{h~^) for the number of resonances close to non-critical energy levels in 
domains of width h and a Breit-Wigner approximation formula for the derivative of the 
spectral shift function. 
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1. Introduction 

The resonance theory for the Schrodinger equation has been developed following several 
approaches. Among them we can mention the analytic dilation (see [1]) or the analytic distor- 
tion (see [21]) and, in the semi-classical regime, that related to the work of Helffer-Sjostrand 
[20]. In [18] Helffer-Martinez showed that the different definitions give the same resonances 
when one can simultaneously apply them to an operator. For the three dimensional Dirac 
operator, Seba [41] defined the resonances as complex eigenvalues of the operator obtain by 
a complex dilation. Applying the approach of Helffer-Sjostrand [20], Parisse [29] has studied 
the Dirac resonances in the semi-classical regime, with some scaling functions. The last two 
works deal with analytic perturbations near the real axis. 

The concept of the spectral shift function has been introduced by Lifshits [25] in connection 
with the problems in quantum statistics and solid physics. Thereafter, a mathematical theory 
of the spectral shift function has been constructed by Krein [24]. Moreover, in [3] Birman- 
Krein found a connection between the scattering theory and the theory of the spectral shift 
function. A detailed presentation of the theory of the spectral shift function can be found 
in [44]. For a survey concerning the spectral shift function (SSF) for Schrodinger and Dirac 
operator or the asymptotic expansion of this function, we refer to Robert [35] and to the 
references given there. 
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A representation of the derivative of the scattering phase in terms of the resonances has 
been estabhshed for the Schrodinger operators. Such representations have been successively 
obtained by Melrose [26] for obstacle problems in the high energy case, by Petkov-Zworski 
[31], [32] for "black box" scattering with compact perturbations in the classical and the semi- 
classical cases and by Bruneau-Petkov [8] for long-range perturbations in the semi-classical 
"black box" framework. The results in [8] have been generalized by Dimassi-Petkov [13] for 
non-semi-bounded Schrodinger type operators. As a by-product, they prove a Weyl type 
asymptotic for the the scattering phase. Moreover, Weyl asymptotic can also be obtained by 
representation of the derivative of the spectral shift function involving the trace of the cut-off 
resolvent (see Robert [36], Bruneau-Petkov [7] and Nakamura [28]). 

Concerning the Breit-Wigner approximation for the derivative of the spectral shift function 
in the Schrodinger case, same kind of results have been obtained in a particular semi-classical 
set-up by C.Gcrard-Martinez-Robert [16] for short range potentials on M" and by Petkov- 
Zworski [32] for a general compactly supported perturbation (see also [7], [6]). 

For Dirac operators, Bruneau-Robert [10] established an asymptotic expansion of the scat- 
tering phase ,s(A) and their derivatives in the high energy regime and in the semi-classical 
regime for A in a non-trapping energy interval. For an interval / c] — mc^,mc^[ with non 
critical extremities, Helffer and Robert in [19] gave an asymptotics of the number of the eigen- 
values in / for scalar potentials. Nevertheless, we are not aware of works which deal with the 
link between the derivative of the SSF and resonances for the semi-classical Dirac operators 
(in the spirit of Petkov-Zworski [32] and Bruneau-Petkov [8]), neither of papers giving the 
Weyl asymptotic of the spectral shift function for Dirac operators in any interval /. 

The purpose of this work is to extend the definition of resonance for analytic perturbations 
outside a compact set. We define the resonances for the semi-classical Dirac operator as the 
discrete eigenvalues of the non-selfadjoint operator obtained from the Dirac operator H by a 
general class of complex distortions of M^. We prove that the resonances are independent of 
the distortion (see Section 4). We establish an upper bound for the number of resonances in a 
compact domain Q (see Section 5). The second goal of this work is to obtain a meromorphic 
continuation of the derivative of the spectral shift function ^(A, h) related to the resonances for 
the semi-classical Dirac operator (see Section 6). This last is closely related to trace formulae 
(see [8], [38], [39], [31], [32], [40]) and to resonance expansions (see [42], [11]). Thereafter, in the 
case where the potential is an electro-magnetic potential, we deduce a Weyl-type asymptotic 
of the spectral shift function (see Section 7). As a by-product, we obtain an upper bound 
0(/i^^) for the number of resonances close to non-critical energy levels in domains of width 
h (see Subsection 8.1), as well as a Breit-Wigner approximation for the derivative of the SSF 
(see Subsection 8.2). 



with domain D{Ho) = H^{R^) (g) C H = L'^{R^) (g) C^, where /i \ is the semi-classical 
parameter, m > is the mass of the Dirac particle and c is the speed of light. The quantities 



2. STATEMENT OF THE RESULTS 



We consider the selfadjoint Dirac operator 
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Oil, a2, 0:3 and (3 are 4x4 Dirac matrices satisfying the anti-commutation relations 

, . ( aiaj + ajtti = 2Sijl4, for i,j = 1,2,3, 

^ ^ \ aiP + Pai = 0, for z = 1,2,3, 

and = I/i. Here In is the n x n identity matrix. For example, we choose the standard (or 
Dirac-Pauli) representation of these matrices 

f ai\ „ f h 



where {(^j)i<j<3 are the 2x2 Pauli matrices: 

01\ f -i \ f 1 



^Ti = I ^ Q J , C72 = 1^ . ; ' - V -1 

Remark 1. Most calculations with Dirac matrices can be done without referring to a particular 
representation (see Appendix l.A [43, Chap. 1]/ 

Let Hi = H -.= 11(1 + V , where V is the multiplication operator by a 4x4-matrix potential 
V . We suppose that V G C°°(R^) and satisfies the following assumption 



(Av) : V is Hermitian on and has an analytic extension in the sector 

(3) Ce,o := {2 G C^|Im(z)| < e|Re(z)|, \Re{z)\>Ro}, /or < e < 1. 
Moreover, for x £ C^^ it satisfies 

(4) \\Vix)\\=Oi{x)-'), 5>0, (x) = (l + |x|2)i. 

The free Dirac operator Hq has only essential spectrum dessiHo) =] — 00, —mc^] U [mc^, +oo[. 
Under the assumption (Ay) the operator Hi is a selfadjoint operator. Using Weyl theorem, 
we have aess{Hi) = aessiHo)- 

For eeDe:={ee C, \e\ < re := ;7=f }, we denote 

Hi^g = He := UeHoUg^ + UeVU^^ = Ho,9 + UgVUg\ 

where Ug is the one-parameter family of distortions defined below (see Section 3). 
For 6*0 fixed in D+ := A n {6* G C, Im(6') > 0}, we define 



re,:={±cJj^^^ + m^c^eC, Ag[0,+oo[}, 
and ' 

Seo-={z^ [j^e; arg(l + 0) < arg(l + ^o), ^^^^ < m^^^ }• 
6>eD+ I I I 01 

The square root ^/z is defined such that for z G C\] — oo,0], Re(y^) > 0. 

For 9 G -D^, arg(l + ^o) < arg(l -|- 9), < ]T+^' prove that the spectrum of Hg is 

discrete in Sq and independent of 9 in Sg^ . This justifies the following definition 

Definition 1. The resonances of H in SggUM are the eigenvalues of Hg^^. The multiplicity of 
a resonance z is the geometric multiplicity of z considered as an eigenvalue of Hg^ . We will 
denote Res(iJ) the set of resonances. 
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The most important advantage of this definition is that the resonances can be computed 
by solving a non-selfadjoint eigenvalue problem. 

Remark 2. The resonances of H in {z G C; Re(2;) g] — m(P,m(P[} are the real eigenvalues 
ofH. 




Fig. The set Se^ 



Now, we would like to find symmetry properties which makes it possible to limit our study 
of the resonances in a domain fl which satisfies (A^j), with 

(A^): ft is an open simply connected and relatively compact subset of {z G C;±Re(2;) > 
mc^} such that Q, PI {±Im(2;) > 0} 7^ and there exists Oq G such that fl n = 0. 

Proposition 1. Let H" be the self adjoint Hamiltonian 

H- =Ho-UcV{x)U-\ 

■where Uc = ij3a2 is an unitary matrix 4x4 and V is the conjugate of V. Then the following 
assertions are equivalents: 

(i) The complex value z is a resonance of H. 

(a) The symmetric of the conjugate —z is a complex eigenvalue of UgH~U^^ . 

Moreover, the multiplicity of z is equal to the multiplicity of —z considered as an eigenvalue 
ofU-eH-m\ 
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Proposition 2. Let H be the self adjoint Hamiltonian 

3 

H = -ic^a'jd^j + P'mc^ + V, 
i=i 

where a'l = —ai, = 012, a'^ = —03, = P are matrices which satisfy the anti- commutation 
relations (2) andV is the conjugate ofV. Then the following assertions are equivalents: 

(i) The complex value z is a resonance of H. 
(a) The conjugate z is a complex eigenvalue of UgHU^^ . 

Moreover, the multiplicity of z is equal to the multiplicity of z considered as an eigenvalue of 

Using the same type of approach as in [39], we construct an operator Hj Q for j = 0,1, so 
that 

Hj^0 - Hjfi = K.J = 0(1), has finite rank 0{h~^), 
and \\{Hjfi — 2:)~"^||=0(1), uniformly for z G (see Subsection 5.1). 

Using this construction we establish an upper bound of the number of resonances: 

Theorem 1. (Upper bound) Assume that V satisfies the assumption (Av) with S > 0. Let 
Q be a complex domain satisfying the assumption (A^), then 

#Res{H) nn< c{9.)h-^. 

For a pair of self-adjoint operators {Hq,Hq + V) where V satisfies the assumption (Av) 
with (5 > 3, (see [10], [36], [35]), the spectral shift function i{\,h) is a distribution in P'(]R) 
such that its derivative satisfies: 

(5) (^'(A,/i),/(A))^,(K),i,(M) = tr(/(i7i) - f{Ho)), /(A) G C^iR)- 

By the Birman-Krein theory, the SSF is in L^^^^(M) and coincides with the scattering phase 

(see]44]). 

Our principal result is a meromorphic continuation of the derivative of the spectral shift 
function ^(A, h). 

Theorem 2. (Representation formula) Assume that V satisfi.es the assumption (Av) with 
(5 > 3. Let be a complex domain satisfying the assumption (A^) and W ^ O, be an open 
simply connected set which is symmetric with respect to the real axis. Assume that 7 = W^flR 
is an interval. Then for all \ E I we have the representation: 

(6) e{Kh) = -lmri\,h)+ + E '^«'(^)' 

TT ^ — ' TT A — wr — ' 

«;eRes(iii)nQ «)eRes(i?i)n/ 

ImtiJ^O 

where r{z,K) = g{z,h) — g{z,h), g{z,h) is a holomorphic function in O which satisfies the 
following estimate: 

(7) \giz,h)\<C{W)h-^ zeW, 

with C(W) > independent of h €]0, ^o]- Here Syj{-) is the Dirac mass at w eM. 



6 ABDALLAH KHOCHMAN 

Remark 3. This theorem can be extended to the pair operators {Hq + Vi, Hq + V2), where 
Vi, V2 are two 4x4 Hermitian potential matrices satisfying (Ay) with 6 > and V = V2 — V1 
satisfies the assumption (Av) with 6 > 3 (see Theorem 5). 

As a corollary of the last theorem, we have a Sjostrand type local trace formula (see Theo- 
rem 6). 

Now we discuss a Weyl type asymptotic of the spectral shift function ^(A, /i) in the case 
where V is an electro-magnetic potential 

3 

(8) V{x) = e(-a ■ A + v){x) = - ^aj ■ eAj{x) + ev{x), 

satisfying the assumption (Ay) with 6 > 3. Here e < is the charge of the Dirac particle. 
We assume that, the electric potential ^(a;) = ^ ^+(|^)-^2 ^ ^0^^ ^ where 1;+, v- are C°° 
scalar functions satisfying 

(9) \e{v+ - v^){x)\ < 2mc^, 

and A = (Ai, A2, A3) is a magnetic vector potential where Ai, A2, A3 are C°° scalar 
functions. 

For any (x,^) G M^, the semi-classical symbols of H^, u = 0,1 are the matrices 

(10) V,y(x,(,) = a ■ {c£, — i'eA{x)) + Pmc^ + i'ev{x), a = (ai, 02, as), 
which is Hermitian and has two eigenvalues 

(11) H^ix,0 = ± (\c^ - veA{x)f + (mc2 + vU,v+ - v^)f^ + vU,v^ + 

of multiplicity two. The function B.^{x, ^) is the Hamiltonian for a relativistic classical particle 
and the other fl^{x,^ can be considered as the one for its anti-particle (see [14], [43], [46]). 
Moreover, from (9) the two Hamiltonians H^{x,$,) are smooth functions. 
For v = 0,1, the matrix 

is the orthogonal projection on the eigenspace £^{x,^) of Pi/(x,^) corresponding to the eigen- 
value H^{x, ^). 

Definition 2. A real X is a noncritical energy level for Hi if for all (x,^) G M®, with 
H^{x,0 = K we have V^,^i/f (x,C) ¥= 0. 

Theorem 3. (Weyl formula) Assume that the potential V is an electro-magnetic potential 
given by (8) and satisfying the assumption (Av) with S > 3. For all X, Xi noncritical energy 
levels for Hi such that ^mc? ]Ai, A[ and h g]0, /io[, we have the asymptotic expansion 

(13) e(A, h) - e(Ai, h) = w{X, Xi)h-^ + 0{h-^). 

Here the 0{h~^) is uniform for X (resp.Xi) in 0. small interval I2 (resp. Ii). The first term 
w{X,Xi) G C°°(/2 X Ii) is given by 

u;(A, Ai) = ■u;(A) - ■u;(Ai), 
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with: 

(14) w{\) = ^ f W^{\v+,V-)-W+{\{),^)-W-{\v+,V-) + W-{\,^,{)) dx 

3 

where W±{X, a, h) = {{x - - (mc^ + ^) ' . 

The Theorem 3 can be extended to the pair operators (i?o + ^ij -^0 + ^) (see Remark 10). 

Remark 4. The two formulae (13) and (14) give in particular a Weyl type asymptotic of 
the counting function of the number of eigenvalues of Hi between two values in the interval 
] — m(?,m(?[ . In the case of scalar potential v {vj^ = V-), this result was proved by Helffer- 
Robert [19] without the analytic assumption at infinity. 

To prove Theorem 3 we construct, in Appendix A (see Theorem 9), a parametrix at small 
time of the propagator of the Dirac equation in an external electro-magnetic field (see also 
Yajima [46] for a scalar electric potential cases). 

As a direct result of the last theorem we deduce an upper bound 0{h~^) for the number 
of resonances close to non-critical energy levels in domains of width h (see Proposition 11) 
and a Breit-Wigner approximation for the derivative of the spectral shift function ^(A, h) (see 
Theorem 8). 

3. Distortion for the free Dirac operator 

In this section, we start with the definition of the deformation for the free Dirac operator by 
analytic distortion (in the spirit of Hunziker [21]) and we will calculate the essential spectrum 
for the distorded free Dirac operator. Here, h does not play any role, and can be taken equal 
to 1. Let us now introduce the one-parameter family of unitary distortion 



Uef{x) = Jl^^^fmx)), /G(5(M3))^ 



where (j)0{x) = x -\- 6g{x) and g : i — > is a smooth function. Let J^g(2.)=det(I-|-6'V(7(x)) 
be the Jacobian of (t>e{x). 

We suppose that g satisfies the assumption 

' (i) sup^gK3 \\^9{x)\\ = < +00. 

(Ag) < (ii) g{x) = 0, in the compact B{0, Rq), (see (3)). 

(iii) g{x) = x, outside of a compact set K{d B{0, Rq)). 
Lemma 1. For 9 e] — M, M[, 17$ can be extended as an unitary operator on H. 
Proof. Since |6'| < M, we have ||0V£?(a;)|| < 1, and 

oo 

ivMx))-' = ii+evg{x))-' = Y,i-mon^9{x)r. 

n=0 
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The function (f)0{x) is injective and ^^(R^) = R^, consequently ^e(x) is a diffeomorphism from 
R^ to R^. The inverse of Uq is given by 

The lemma follows from the relations 

UeU^^ = U^^Ue = /(l2(k3))4 and ||C/,/||w = \\U^'f\\H = \\f\\, V/ G iL\R^))\ 

□ 

Definition 3. We denote by B, the space of functions f = (/i)i<i<4 such that fi{x) has an 
analytic continuation in Cefl and lim \z\''fi{z) = 0, for all k EN and e g]0, 1[ (see (3)). 

Lemma 2. The subspace B is dense in 7i. 

Proof. The subspace B contains vectors of Hermite functions and the linear combinations 
of Hermite functions are dense in L^fR^). □ 



Proposition 3. Let be D^^m = -De H {0 G C; \9\ < M}. We have the two assertions: 
(i) For all f E B, 9 E D^^m ' — Ugf is analytic, 
(a) For all 9 G D^^m, UqB is dense in Ti. 

Proof. In the order to prove (i), we show that 9 i — >■ {Uof,g) is analytic for all g E H. Let 
i? > 1 be such that K C B{0,R) = {x, \x\ < R}. 

1 

• In {x, \x\ < R}: since f E B, then 6 i — > J^^^^j^ J^^f {(j)Q{x))g{x)dx , is clearly analytic 
for all g E Ti. 

• In {x, \x\ > R}: we have g{x) = x, consequently </>6)(x) = (1 + 9)x. We remark that 

If l^^l <re = :7=5, < e < 1, then 
According to the definition of B we have 



then, 9 i — Jj^|>^ J^J{^e)g{x)dx is analytic, 
(ii) Let h{x) E (C^(M3))^. We denote 

hk{x) = (^)t j e-''^^-y-'^^y^^'h{y)J^^^y^dy, 

which is clearly in B. 

Using (1)1 /e-*^(^-f-^5{s/))V^^(^)(iy = f e-'^^' dz = 1, we get 

h{x) - hkiMx)) = (^)^ / e-''(M-^-My))\h{x) - h{y))J^^^y)dy. 
The last term tends to when k — +oo. Consequently, we have 

hk o 4>e{x) h{x) , in 7^. 
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□ 

Remark 5. One can always choose g satisfying the assumption (Ag) with M > = :^7f^^- 
In that case, we have -De,M = D^. 

Lemma 3. For 9 € D^, we have 

1 ^ d 

j=i 

where Qe{x,dxj) = ^\a\<i'^oi{x,9)d'^. is such that: 

(i) 9 I — > aaix,9) is an analytic function bounded by 0{9). 

(n) x^a,,{x,9) € (C^iR'^))^. 

In particular 9 i — > Hq^q is an analytic family of type A of domain D(Hq) (see Kato [23], for 
the definition of an analytic family of type A). 



Proof. We denote dj = and we calculate the term UedjU^ ^. 
UedjU^^fix) = Ugdj (jffi(f>oHx)) 



1 ^ 

-2'^^e^3J<t>e{.x)f{x) + ^dkf{x) (s^^^^fc) (</)(?(x)), 

U — 1 



with ^/(x) = (polix), 0ej(x)y 

We remark that (f>g^{x) = outside of the compact set K, then 

3 ^ 

{(j)0{x)) = Y^djf{x), outside of K. 

k=l 

Let xk £ C'o°(^^)' — Xk < 1, be equal to 1 on if and outside a compact set which 
contains K. We have 

-, , 3 

UgdjU^'fix) = --J-;a,J^^(,)/(x) + — -5,/(x)(l-XK) + E9fc/(x) (d,4>e'k) iMx))XK. 

k=l 

Since djJ^^^^x) has a compact support, then 

(15) UedjU^^ = j^^dj + qeix, d,.) 

with qe{x,dxj) satisfying the hypothesis of Qff{x,dxj). 

Now we just have to multiply (15) by —icaj, to sum on all values of j and add /3mc^ to 
both hands. The estimate aa{x,9) = 0{9) is clear using that Qo{x,dx^) = and the analytic 
property of 9 i — >■ aa{--,9). □ 
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Lemma 4. Let Pg = TTo{—ic^^ ^jlir) + Pmc? . Then 



a{Pe) = aessiPe) = rg = {zeC-z = ±c{ ^ + m^c^)^ , a € [0, +oo[}. 
Proof. Let be the Fourier transform and 

j 

mc^h (i^)((^i6 + cr26 + c^36) 

where ^ = (^1,^2)^3) € and aj^j is the multiplication operator by the 4x4 matrix Oij^j. 
Here (7i,a2,crs are the 2x2 Pauli matrices. The spectrum of Pg coincide with the spectrum 
of the multiplication operator K{6). We easily prove that 

and we deduce the lemma. □ 
The principal branch of the square root function is holomorphic on the set C\] — oo,0]. Let 
= {z= jy:^ + m^c^, 9 e D^, \ e [0, +oo[}. Since, 

Sd, C]0,+oo[el"5'5[, 
the square root z 1 — ^ z^ is holomorphic on Sd^- 

Lemma 5. For Hq q, Pq, defined as above, we have aess{HQfi) = (Tess{Pe)- 

Proof. We want to use Kato's theorem [23, Th.4.5.35]. For A > 1, A G M and Qe defined 
in Lemma 3, we have 

{H^fi - i\) = (1 + Qe{Pe - iX)-^){Pe - iA). 

Since {Pe - iX)'^ G C{n, (Hy) and Qg{Pe - iX)'^ = 0(f ), we obtain that iX G p{Ho^e) = 
C \ u{HQfi). To apply the Kato Theorem, it is enough to show that 

(16) {Hofi — iX)~^ — {Pq — iX)~^ is compact. 

Using the resolvent equation, we have 

{Ho,e - iX)-^ - {Pe - iX)-^ = {Ho,e - iX)-^Qe{P6 - iX)-\ 
with Qo{x^ dxj) = ^\a\<i 0'a{x, d)dx- compactly supported. Since the operaror {Ho^Q—iX)''^Qe 
is bounded and lsupp(Qe)(-Pe ~ iX)~^ is compact, the assertion (16) holds. □ 

4. Definition of resonances 

In this section we distort the perturbed Dirac operator H = Hq + V, where the potential 
V satisfies the assumption (Ay) and we define the resonances for the semi-classical Dirac 
operator. 

The distorted Dirac operator is denoted by 

He = UeHoUg' + UeVU^^ = Ho,e + V{Mx))- 
Proposition 4. We suppose that the potential V satisfies the assumption (A-y) , then 
(i) E I — ^ He = -ffo,6» + V{<p9{x)) is an analytic family of type A. 



a{K{9)) = aess{K{9)) = Fg = {z G C; z = ±c( ,^ , + mV)2 , A G [0, +oo[}, 
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(a) aess{He) = Fe. 

Proof. The assertion (i) is clear since Hq^ is an analytic family of type A and V satisfies 
the assumption (Av)- 

Now, we prove (ii) as in the proof of Lemma 5. For A 3> 1, iA € p{Hs) and 

(17) {He - - iHo,e - iXy' = {He - iX)-W{Mx)){Ho,e - ^A)-^ 

Since the operator V{(j)e{x)){H()^o — is compact (see (Ay)) and the resolvent {Hg — iX)~^ 
is bounded, the difference {He~iX)^^ — {Ho^e — iX)~^ is compact. According to Kato's theorem 
[23, Theorem. 4. 5. 35] and to the lemmas 4,5, we obtain (ii). □ 
We denote 

E = {z G C; lm{z) > 0, Re{z) > -mc^} U {z e C; Im {z) < 0, Re{z) < mc^} \ a{H). 

Theorem 4. With the notations used above and taking Oq G -D+ = D {Im(^) > 0}, we 
have: 

(i) For all f,g G B, the function: z G E i — ^ Mf^g{z) = {{z — H)~^f,g) has a meromorphic 

extension on Sq^^. 
(ii) The poles of Mf^g{z) are the eigenvalues of He^. 
(Hi) These poles are independent of the family Ue^- 

(iv) aci{Hef^) n S = 0, where ad{Hef^) is the discrete spectrum of the operator He^. 
Proof, (i) Since Ue is unitary for G M, 

Mf^g{z) = {{z-H)-'f,g) = {{z-Her^Uef,Ueg). 

We denote 

(18) Mf,g,e{z) = {{z-He)-'U0f,Ueg), for 9 e D,. 

According to (i) of Proposition 4 and to the definition of Ue, the functions i-^ {z — He)^^ , 

6 ^ Ugf and 6 i— > {ip, Ugg) are analytic on for all tp G Ti and any z G S. 

Thus, for z G S, the function 6 i-^ Mf^g^g{z) is analytic on D^. Since Mf^gfi{z) is independent 
of on the real axis and according to the uniqueness of the prolongation, its extension is 
independent of 6. 

Now, we fix 6q G Df. Since Se^ n Uess{Heo) = 0, the function z G S i — ^ Mf^gfi^{z) has a 
meromorphic extension in Se^- 

(ii) First, let z G Se^ be a pole of Mf^g{z) which is equal to Mf^gfi^{z) for G Df. Then 
z G (Td{He,^ n S'e,, (see proof of (i)). 

Now, let e G (Jd{HeQ) n Sg^. There exists u £ 7i such that ||tt|| = 1 and Hg^u = eu. Let 

7 be a small disk centered at e such that jD a{He^) = {e} and F be the positively oriented 
boundary of 7. 

Let us introduce the projector 

n = -^ [ {z - He,)-^dz; Uu = u. 

Since Ug^B = H = Uq^ B (see Proposition 3), there exist /„, gn £ B such that 

\u - Ueofnl < - and \u - UEgn\ <-, n G N. 
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Therefore, as n goes to infinity, we have 

^Jii^ - Heo)~^Ueofn, U0-^gn)dz = ^Jii^ - Heo)~^u, u)dz + o(l) 

= {Uu,u) + o{l) 

= \Wf + o{l) 

= l + o(l), 

and then, 



^- J^{{z - H)-^U,9n)dz = 1 + o(l). 



2i7r 

So that Mf^^g^{z) admits e as a pole in 7. 

The assertion (iii) results from (ii) because Mf^g{z) is independent of Uq. 

(iv) If there exists z G aaiHg^) n S , then z is a pole of {{z — H)~^f,g), for f,g & B, but 
((z — H)~^f,g) is analytic on this domain. We conclude that such z does not exist. □ 

Remark 6. (i) It results from (ii) of Theorem 4 that for all 9 G -D^, the discrete spec- 
trum ad{H) is a subset of ad{Hg). 

(ii) The previous theorem justifies the definition of the resonances (Definition 1) and using 
Lemma 4, Ho has no resonances. 

Remark 7. If 6 G D^, then its conjugate 9 G D^. Repeating arguments of the proof of Theorem 
4, we have 

(i) The function 9 1 — >■ Mj^g Q^z) has a analytic extension for 9 E D^. 

(ii) The function 2; G S 1 — ^ M^^ q^z) has a meromorphic extension on Sg^, 
where 

t = {zeC; lm{z) > 0, Re{z) < mc^} U {2; G C; lm{z) < 0, Re{z) > -mc^} \ a{H), 
and Sg^ is the symmetric of Sg^ with respect to the real axis 

Seo = {^^ U^e; arg{l + 9)<arg{l + 9o), < }• 

eeD+ I I I 01 

Consequently, we obtain (see Theorem 4): 

1) The poles of Mf^g{z) in Sg^ are the eigenvalues of Hg^. 

2) These poles are independent of the family Ug^ . 

3) ad(i/ejns = 0. 

The assertions 3) and (ii) prove that the operator Hg has only discrete spectrum in Sg^ . 
Proof of the Proposition 1. We consider the anti-linear application on H, 

C : I — > Ucip = i(3a2'(p- 

Then, we have 

CHgC-^ = -Hf^g + CVo(f>g{x)C-^ 

= -(H^^g-UcVo<j>g{x)U-'). 
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Using that V is analytic, we get V o (^^{x) = Vo <^§{x). Then, 

CHeC-' = -{Ho^s-U,Vo<Pg{x)U-') 

= -{H,^e-UcUeV{x)U^'U-'). 

1 

We recall that U0f{x) = J^^^^^f{(f)o{x)). Since UcUg = UqUc, we obtain 

CHeC-^ = -[H,^g-UgUcV{x)U-%') 

= -Ug {Ho - UcV{x)U-^) Ug^ = -Hg. 
Consequently, C {Hq — z) = — {H^ + z), and the property follows. □ 

Proof of the Proposition 2. 

By definition of Hq, we have 

(19) Hi^ = Ug {Ho + V- z) 

Using that ai = ai, a2 = —0:2, 03 = 0:3, /3 = P, we find 

^9 (9 

Hq = ^c^^Oj— h Pmc^ = ~^c^^a^ — 1- (3'm(? 

j=l j=l 

and 

Ho + V = H. 

Using the last relation and the equation (19), we obtain Proposition 2. □ 

Finally, the study of resonances in a domain of the complex plan C is reduced to study 
the resonances in Qr\{z G C, Im(^) < 0}, with Q, satisfying the assumption (Aj!j) (see Fig.l). 

5. Upper bound for the number of resonances 

In this section, we establish an upper bound of the number of resonances in a compact 
domain Q. For this purpose we construct an operator Hg : D{Hq) Ti. with some properties 
(see Proposition 6). According to the Section 4, it is sufficient to treat the case where Q, 
satisfies the assumption (A^). 

We shall use the theory of ^-pseudo-differential operator (see [12], [34]). Let m he an order 
function on M^" (i.e. there are Co, A^o > 0, such that m,{x) <Co{x — y)^^m{y)). The space 
SP{m) is the set oia{x,^;h) G C'=^{R^'^)®& such that for every a G N^'', there exists > 0, 
such that 

\\Vl^a{x,i-h)\\<Cani{x,Oh-^. 
For a symbol a{x,^;h), we define the Weyl quantization, a^{x,h'Vx',h) := Op^(a) by 

where u{x) is in the Schwarz space. 
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5.1. Construction of Hq. We follow the approach of Sjostrand [39]. Let be a complex 

domain satisfying the assumption (Aj'j) and ij) € Co°(M^) be such that 'ip{x) > 0, ip{x) = 1 if 
\x\ < 1 and ip{x) = if |x| > 2. We recall the notations of Section 3: </>6)(x) = x + 9g{x) with 
g{x) = in the compact set -6(0, Rq) C K and g{x) = x outside K C -8(0, ao) where ao > 
is sufficiently large. 

Using Lemma 3, the semi-classical principal symbol of Hg is given by: 

1 ^ 

heix,^) = a ■ Ce{x,C) + mc^P + V{(l)g{x)) - c-^aj • J^g(^)djJ^g{x), 

with 



0(a;,0 = (Ce,i(a^'0> C.e,2{x,0^ C(?,3(a;,6) and Ce,i(a;,0 = c^^k \dj(t)QjA {(^eix))- 



k=l 

„2, 



For all (x,^), the matrix M = a ■ Ce{x,0 + tic /?, has two eigenvalues 

A± = ±VCe(x,e)2+mV. 
Consequently, there exists an invertible matrix U such that 

where 



U 



^2 



with cr = (ci, (72, era) and ('7j)i<j<3 the 2x2 Pauli matrices. 

One can easily prove that the norms of U, U"^ and their derivatives are bounded in the 
following way: 

\\d^cf^uix,o\\ < c{0-^ 

(20) \\d^d^U-\x,m < CiO'^ Va,/3eN. 
Applying i7~^(resp. U) on the left (resp. on the right) of ho{x,^), we obtain 

U-^he{x,OU = de + Ve{x,0, 

where Vo{x,S,) = JJ-'^ (v{4>g{x)) - c\ Y^^=x «i ' '^4^g(x)^j'^M^)) ^- ^^^^ ^'^^ 

is compactly supported _7 = 1,2,3, (see Lemma 3) and the matrix V{(j)e{x)), U, U~^, and their 

derivatives are uniformly bounded, then V$(x,^) and their derivatives are uniformly bounded. 

In order to construct Hg, we introduce an intermediate function f{x,^): 

We denote the diameter of J7. Let us choose Po > and Co > 0, sufficiently large 
such that 

(21) V^GM^; sup\\Ve{x,0\\ + M <l\>^^-iCo^P{-l-)\ 

xe MS ^ Po 



y(Re(A±))'+ (lm(A±)- Co^{£ ' 
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We prove (21) considering the two cases: 

• For 1^1 > (3q, with /?o > 0, sufficiently large we have 

sup \\Veix,m + M<hM^^)\- 

• For 1^1 < /3o, since is bounded, we choose Co > 0, sufficiently large such that 



1 



sup ||F,(x,OII + M<^|Im(A±)-Co|. 

For |x| > cuo > 0, sufficiently large we have (^^{x,^) = and = icy^ (i+e)^ m^c^. 
Since the domain CI satisfies the assumption (A^^), we have 

min{dist(n, A^), dist(n, A^)} 7^ 0, 
hence we can choose ao > sufficiently large such that 

(22) V|a;| > ao, \\Ve{x,0\\ < \ dist(n,A^) := ^ min{dist(n, A+), dist(n,Ag )}. 
Now, we define f{x,$,) in the following way: 

(23) f{x,o = coi;{^mh- 

Lemma 6. The matrix he{x,^) — if{x,$,) — z is invertible for all z E and satisfies 

(24) \\d^d^^iheix,0-if{x,0-^r'\\<C{0~'-^, Va, /3 G N. 

Proof. Applying C/~^(resp. U) on the left (resp. on the right) of hQ{x,^) — if{x,^) — z, 
we obtain 

{he{x, - if{x, O-z)U = d0- if{x, + Veix, 0- 
1) Let us prove that the symbol a := dg — if{x,^) — z + Vo{x,$,) is invertible. 
• For l^l < ao, 

a = {de-iCo^P{j-^))(^h+(^dg-iCo^{^^)j {Ve{x,0-z)y 
According to (21), we have 



(^d,-iCoV(^)) (V0{x,o-z)\\<^, 



thus a is invertible and satisfies 



(25) \\a-'\\ < 2|| (^de - ^CoV'(^)) II < CiO''. 

• For \x\ > ao, we have A^ = icy^ (i+ep' ^^c^- Since f{x,^) > 0, we have 
|A+ -iz + if{x, 0)\ > dist(n, A±) > C(0 > 0, 

and 

|Re (Ag -{z + if{x, 0)) I = |Re (Ag - z) \ > dist(ll, A^) > C(0 > 0. 
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Since 

a = {de- if{x, - z) (/4 + ide - ifix, - z)-^ Ve{x, o) , 

and 

II {de - if{x,i) - z)-^ Ve{x,m < ^, (see (22)), 
the matrix a is invertible and 

(26) \\a-^\\<2\\{de-if{x,i)-z)-^\\<C{i)-\ 

2) According to 1), the matrix J7^^(/i5)(x, ^) — i/(x,^) — z)U is invertible. From (25), (26) and 
(20), we deduce that the matrix hg{x,^) — if{x,$,) — z is invertible and 

Uh0{x,O-if{x,O-zr^\\ = \\U-' (de-if{x,Oh-z + Ve{x,Oy^U\\ 

< ||C/||||C/-i|||| (d0-if{x,Oh-z + Veix,Oy^\\ 

(27) < C{0-'. 

This gives (24) for a = /3 = 0. Using (20) and (27) we obtain (24) for {a,P) G by induction. 

□ 

We denote Hg = Hg + T, with T = Op%{-if{x,(,)), where f{x,£,) is defined in (23). It's clear 
that the semi-classical principal symbol of {Hg — z) is 

Proposition 5. If h > is small enough, the operator {z — Hg) is invertible for every z G O 
and, for every N G'N its inverse satisfies: 

{z - Hg)-' = On{1) ■■ D{H'') ^ D{H''+'), 
uniformly for z € fi. Here D{H^) designates the domain of H^ with the convention D{H^) = 

n. 

Proof. Let us prove that the operator {z — Hg) is a Predholm operator of index 0. We 
have 

{z - Hg){z - Ho^g)-' = {z - Ho,g + Ho^g - Hg){z - Ho,g)-' 

= I-{f + V{ct>g{x))){z-Hofi)-\ 

Since the right-hand side is a perturbation of the identity by a compact operator and 

{z-Ho^g)-' : (l2(IR3))^ ^ D{H) is invertible, 
the operator {z — Hg) is Predholm of index 0. Consequently, it is enough to show that 

(28) \H\h'^+^) < C\\{z - Hg)u\\l^H^y for u G D{H''+'). 
According to Lemma 6, the symbol qo = is well defined and satisfies 

Moreover, having 

W.dl{a^^)\\<C{C,^'-^ 
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the composition theorem of /i-pseudo-differential operators imphes 

Opt{qo)Opt{aH,) = Opt{r) 

where (r — 1) is in the space of symbols S^\h). In particular the operator 

Op^(r) : D{H^+^) ^ D{H^+^), ViV € N, 

is invertible for h small enough, then (28) follows. Therefore the operator {Ho — z) is also 
invertible and we have 

{z - He)-' = On{1) : DiH"") ^ D{H''+'). 

□ 

Proposition 6. There exists Hg : D{H) \ — > Tl, with the following properties. 

The difference K := Hg — Hg is of finite rank 0(h~^), has a compact support in the sense that 

K = XiKxi for some xi £ C^(M-^) and 

K = 0(1) : D{H^) I — > D{H^) ViV, M G N. 

Moreover, for every N Ef^, we have 

{He - z)-' = 0(1) : D{H^) ^ D{H^+'), 

uniformly for z E Q. 

Proof. (We again use all the previous notations) We define 

He := He + xi^'Xi = He + xi^Xi - T, 

with xi{x) = V'Cii^) 

T := x(-/i' A + x^)T = x(-/i' A + x')Opl{-if{x, 0) 
where x ^ C'q°(R) is such that: 

x{^ + x^) = 1 on the support of f{x,$) (see (23)). 
By functional calculus (see [12]), we can prove that 

(29) He~He = f- xiTxi = 0{h^) : D{H^) i — > D{H^), V M, N eN. 
The last lemma, the formula (29) and 

{He - z)-' = {He - z)-' (l + {He - He){He - z)-")'^ 

yield for all N e N 

{He - z)-' = 0(1) : D{H^) ^ D{H^+'). 

According to the fact that x(~^^A + x^) is of finite rank 0{h~^), to the fact that Weyl 
quantification Opf^{—if{x,£,)) is bounded and to the definition of xi, the operator 

K:=He-He = XI (x(-^'A + x'')Op%{-if{x, 0)) Xi 
is of finite rank 0{h~^) and compactly supported. □ 
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5.2. Upper bound for the number of resonances. In this section we establish the upper 
bound of the number of resonances given in Theorem 1. 

Lemma 7. Let p> he an open complex relatively compact subset ofC and Hg be define as 
above. There exists g satifying (A.g) such that for h small enough and z G rifl {Im z> p> 0}, 
we have (z - Hg)'^ =0(1). 

Proof. We again use the notations of Section 3: 4'0{x) = x + Og{x) with g{x) = in the 
compact B{Q,Rq), and the notations of Subsection 5.1 concerning U, U~^, dg and 

V^(x,^) which satisfy 

U-^he{x,i)U = de + Ve{x,i). 

The matrice /i6)(x,^) is the semi-classical principal symbol of Hg. 

According to Section 4, the resonances are independent of the family Uq. Then we can 
assume that g{x) = in the ball B{0,Rg) D B{0,Ro), with Rg > 0, sufficiently large such 
that 

VxeM^ \x\>Rg>0, \\Vgix,0\\<^. 

Repeating arguments of Subsection 5.1, we can prove that {dg + V^(x,^) — z) is invertible, 
then (/ig(x,^) — z) is invertible and 

\\d^cf^{he{x,0-zr'\\<C{0-'-^. 

Since we have: 

\\d^d^^{hgix,0-z)\\<C{0+'-^, 

the composition theorem of /i-pseudo-differential operators implies 

Op^{{he{x,0 - z)-^)Opl{he{x,i) -z) = l + Oih), 
where 0{h) corresponds to the norm in >C(L^). □ 

Proof of the Theorem 1. 

Let K{z) = K{z — He)~^ with Hg, K defined in Proposition 6. We remark that 

(/ + K{z)){z -Hg) = {z-Hg)+K = Z- Hg. 

Thus, the resonances z G Res(i/)nr2 repeated with their multiplicities coincide with the zeros 
of the function 

D{z) = det{I + K{z)). 
Since K is bounded and is of finite rank 0{h~^), we have 

\D{z)\ < ell-^^^)!!*-^ < e'^°^~\ for all zeU. 
Using Lemma 7, we get {z — Hg)~^ = 0(1) for Imz > p > and z E CI. Since 
(30) iI + K{z))-^ = {z-Hg){z-Hg)-\ 

then 

||(/ + ^(z))-^|| < Ci, lmz>p>0. 
Writing the operator (/ + K{z))^^ in the form 

iI + Kiz))-'=I-K{z){I + Kiz))-\ 
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we obtain 
which implies 



\det({I + K{z))-^^\< e^2/» \ lmz>p, 



\D{z)\ > Ce-^3fe \z enn{lmz> p}. 
Now, applying Jensen's inequality in a slightly larger domain, we obtain Theorem 1. 

6. Representation of derivative of the spectral shift function 



□ 



In this section we prove our principal result given in Theorem 2 and a generalization (see 
Theorem 5). Moreover, we give a Sjostrand type local trace formula. 

The spectral shift function ^{X,h) (g P'(M)) associated to Hq,Hi is defined (see [10], [36], 
[44]) by 

{^'{X,h)J{X)) = tr(/(i7i) - f{Ho)), f e C^{R). 
In the following, we will use the notations: 

Hi = H, Ki ■.=K = Hifi - Hifi --He - He and [a. ]o = ai - do- 
For an integer m > 3, we define the functions: 



(31) a±(z) = (z2 + l)™tr[(if. 
The (j± satisfy the relation 



\H+i)-"'{z-H. 



Jo 



±lmz > 0. 



±Imz > 0. 



(32) a_{z)=a+{z), Im (z) < 0. 

Proposition 7. For a potential V satisfying the assumption (Av) with S > 3, the function 

r ' - 1 

6 I — > {H.fi — i)~™'{H. + i)^'^{z — H.^q)^^ is holomorphic from Df to the space of trace 
class operators. Moreover, for any 9 G , we have 

(33) a±{z) = {z^ + l)-tr [(iJ,, - lY'^iH.^e + - ^i-fiY^ 
Proof. For G M, the operator 

{H. - iy^'iH + iy^'iz - /^.)"^ 

is unitarly equivalent to the operator 

(i7,,-z)--(F.,, + z)— (z-iJ.,,)-\ 
Using the cyclicity of the trace, we deduce 



(34) a±{z) = {z^ + l)™tr {He - i^^iH-fi + i)-™(^ - H., 



±liRZ > 0, 6* G 



According to the proof of Theorem 4, the resolvent {z — H. e) ^ is analytic for 9 G and 
z G n {Imz > 0}. Then, the function 9 i — > {H.^e - i^^iH.^e + i^^iz - H.^)'^ is also 
analytic on Df. 
Now, we treat the difference 

1 



{z-H,e)-\H.^e-i)-'^{H.,e + i)- 



(3c 



= AiBiCi — AqBoCq 

= Ai5i(Ci - Co) + Ai(Si - Bo)Co 
+ (Ai-Ao)SoCo. 
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Clearly, the terms A. := {z - H.^)''^ for Im z > 0, B. := {H.^ - i)-"* and C. := {H.^ + i)- 

are bounded. 

For any integer m > 3, the term 



(36) 



Si(Ci - Co) = (Bi(Ci - C^){xY{hV,r) {{hV,)-^{x)-') , 



is analytic for 9 G with values in the space of trace class operators. This can be proved 
using functional calculus in the framework of /7-pseudo-differential operators (see [12]): The 
first factor Bi(Ci - Cq){x)^ {hV:,)"" is analytic for 6 G D+, the second factor ((W^)-™(x)-^) 
is in the space of trace class operators and its trace norm is bounded by 0{h~^). Then, the 
left-hand side of the equation (36) is in the space of trace class operators and its trace norm 
is bounded by 0{hr'^'). The same argument can be used for the terms Ai{B-i — Bq) and 
(^1 — Aq)Bq, then their trace norm are bounded by 0{h~'^^ 



Since the function tr 



(z - H., 



1 1 



is analytic with respect to 



9 G and independent of 9 on the real axis, the formula (33) follows. 



□ 



Repeating the construction of Hi^, we can construct an operator Hq q : D{Hq) — > Ti with 
the properties of Hq such that the difference Kq := Hq^q — Hq^ satisfies the properties of K\ 
(see Proposition 6). 

Proposition 8. There exists a function a+{z,h) holomorphic in such that for all 
z G O n {Im(2:) > 0}, we have: 



(37) 



a+{z) = tr\{H.^e - z)-'^K.{H.fi - + a+{z,h), 



\a+iz,h)\ <C{n)h-^, zeCl, 
with C(ri) a constant independent of h. 
Proof. For 2; G O fl {Imz > 0}, we have 

(38) {Ke - z)-^ = (H.,e - z)-^ + {Ke - z)-^K.{H.,e - z)-\ 

From the equations (38) and (33), we deduce: 



a+{z) = {{z -i){z + i)r tr [ [{kg - z)-\Kg - i)-^{Kg + i)'™) 

+ {{z -i){z + i)r tr [ ((i?.,e - z)-^K.{H.,e - z)-\H.fi - i)-^{Ke + O"'") 
= A{z) + B{z). 

Starting with the resolvent equation, we obtain: 

{{z -i){z + i)riH,e - ir'^iH^e + ir^'iKe - z)-' 

m 

= (H,g-z)-'-Y,i^ + i)'-\H,e + i)-' 
k=l 

rn 

- {z + i) ^(z - i)''-HKe + t)-^{Ke - i)- 

k=l 
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Using the last equation, the cydicity of the trace and Propositon 6 we obtain 



B{z) = tr 



■flL 

k=i 

m 

{z + i) ^(z - i)'-\H.,e + iT'^iH.fi - i)-^) 



k=l 



tv[{Ke - zr'K.{Ke - z)-')]\ b{z). 



Since the operator {H. q — z)~^ is bounded and holomorphic in J7 by construction, b{z) is 
holomorphic and bounded by 0{h~^). 
It remains to show that 

Aiz) = {{z-i){z + i)rtv[{H,e-z)-\Ko-i)-"'{Ke + ir'^]l 

= {{z - i){z + i)r tr(li5iCi - lo^oCo), 

is holomorphic and bounded by 0{h~^). 

We recall that the terms A. := {H.,e - z)''^ for z e B. := {H.^e - and C. := 

{H.^o + i)"™" are bounded. Using the assumption (Av) with > 3, we treat the difference 
{AiBiCi — AqBqCq) as (35). The only difference is for the term {Ai — A{))B{). We write 

{M-Aq)Bq = {H,^e-z)-HHo,e-H^,e){Ho,e-z)-HHo,e-ir^, 
with HoM-Hi^e = Ho^e - Hi^e + Ko - Ki. 

Then, modulo a trace class operator uniformly bounded, with trace norm bounded by 0{h~^), 
we have 

{A, - Ao)Bo = {Hifi - z)-^ o {{Ho,e - H^,e){Hi,e - iT"') 

o {{H^^e - ir{Ho,e - z)-\Ho,e - i)"'") • 

The second factor {Hqq — Hi q){Hi q — z)^™ is trace class and its trace is 0{h~^), the first 
and the third factor are bounded. Then, the term {Ai — Aq)Bq is analytic for z G O with 
values in the space of trace class operators and its trace is bounded by 0{h~^) and so is the 
difference {MBiCi - AqBqCo). □ 

Lemma 8. For f G Co°(M), we have 

(39) (e',/) = lim^ J f{X)[a+{X + ie)-a.{X-iE)]dX. 

This limit is taken in the sense of distribution. 

Proof. We follow the proof of [13, Lemma 1]. Let / G C^(M), f{z) G 
almost analytic extension of / and 

gix) = f{x){x' + ir. 



^) be an 



Then 



giH.) = ~J dj{z){z'' + iriz - H.r'Lidz), 
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where L{dz) is the Lebesgue measure on C. Clearly 
f{H.) = {H.-i)-^{H. + i)-^g{H.) 

= ~J ^^hz){z^ + ^)"\H. - ir'\H. + i)-"\z - H.)-'L{dz) 

which implies: 

tr(/(iJi)-/(i7o)) = ~jdj{z){z^ + ir 

(40) X tx^{H.-i)-^{H.+i)-'^{z-H.)-^^L{dz). 

We have a±{z) = 0(/i"3|Im .^nd the derivative dj = 0(|lmz|^) for all G N 
(/ G Cq°{R)), so we write the right-hand side of (40) as 

{C'J) = tv{f{H,)-f{Ho)) 

= — — lim [ /" dzf{z)a+{z + ie)L{dz) + ( dzf{z)a-{z — ie)L{dz] 

According to Proposition 8, the functions (T-\-{z + ie) and (T_(z — ie) are holomorphic in 
{z G O; Imz > 0} and {z G fi; Imz < 0} respectively. Applying the Green formula, we obtain 
the lemma. □ 

Before the proof of Theorem 2, let us give the following proposition: 

Proposition 9. (see [38], [39]) Let F[z,h) he a holomorphic function in an open simply 
connected domain O containing a number N{h) of zeros. We suppose that, 

F{z,h)=0{l)eO^^^''^''\ zGn, 

and for all p > small enough, there exists C > such that for all z ^ Qp := fl O {Imz > p} 
we have 

|F(z,^)| >e-^^W. 

Then for all open simply connected subset Q CI there exists g{-,h) holomorphic in $7 such 
that 

N{h) 

F{z, h)= n - Zj)e3^^^''\ dzg{z, h) = 0{N{h)), z e Q. 

i=i 

Proof of Theorem 2. We follow the argument of Sjostrand ([39]). Let 

K.{z) = K.{z - H,gr\ 
As [39, Equation(4.31)] we have the representation 

-ti(^iKe-z)-^K.iH.^g-z)-^^ = ti (^{l + K.{z)y^ d,K.{z)^ 

= dzlog det (l + K.{z)y 
Prom Subsection 5.2 the resonances are the zeros of the function 

D{z, h) = det (^I + Ki{z)^ = 0(l)e'='^ 
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Since the function det(l + Ko{z)) has no zeros in J7 (see (30) and Remark 6), the term 
^zlog det ^1 + Kq{z)^ is analytic and using the Proposition 9, it is bounded by 0{h~^). 
We recall the Res(i7) be the set of resonances of H and let 

D{z,h) = G{z,h) {z-w), 

w 

where, G{z, h) and its inverse are holomorphic functions in J7. Obviously, 

(41) d,\ogD{z,h)=d,\ogG{z,h)+ J2 

Res(if)nn 

Using Proposition 9, we have 

\d^logG{z,h)\ < C{n)h-^, zen, 
where O CC O is an open simply connected set and C{Q) is independent of h. 

Now, we treat the non-holomorphic term in {a^{\ + ie) — a^(X — ie)) when £ — 0, which 



IS 



V , _ ), forA€/. 

^-^ ^ A + le — w \ — le — w I 



w eRes(H)nn 
If Im(w;) 7^ 0, we have 



-1 / 1 1 \ -lm(w) 

lim ' ' — 



2iTr e^o \\ + ie — w X — ie — w J 7r\X — wf 
while for t/; G M we get 

^ lim ( ) = d(X -w) = (5^(A). 

2zTre^o \X + ie-w X-ie-wJ ^ ' ^ ' 

The second hmit is taken in the sense of distributions. 

Lemma 8 and Proposition 8 show that the function r{z^ h) = g{z, h) — 'g{z, h), with g{z, h) = 
a-\-{z, h) + dz\ogG{z, h) + S^log det ^1 + Kq{z)^ is a holomorphic function in Q. and satisfies 
the following estimate: 

(42) \g{z,h)\<C{n)h-^, zeW, 

with C{fl) > independent of h. □ 

Theorem 2 can be extended to a more general situation: 

Theorem 5. Assume that Hi = Hq + Vi, H2 = Hq + V2. The potentials Vi,V2 (resp. 
V = Vi — V2) satisfy the assumption (Av) with S > (resp. S > 3). Let fl be a complex 
domain satisfying the assumption (A^), W <^ Q be an open simply connected and relatively 
compact set which is symmetric with respect to M. Assume that I = W CiR is an interval. 
Then for all X E I we have a representation of the derivative of the spectral shift function 
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associated to the pair {H2, Hi) of the form: 

(43) ei\h) = llmriX,h)+\ J2 Jx'^Z2 + E 

TT L ^ — ■' 7r\A — wr ^ — ' J2 

w&ResiH.)nn weR^s{H.)nI 

where r{z,h) = g{z,h) — g{z,h), g{z,h) is a holomorphic function in O which satisfies the 
following estimate: 

(44) \g{z,h)\ <C{W)h-^ zeW, 

with C{W) > independent of h. Here dw{-) is the Dirac mass at w E^. 

Proof. We denote if2,6» = U0H2Ug^ {Ug defined in Section 3). As in Subsection 5.1, one 
constructs H2fi ■ D{H) 1 — > H with the following properties: 

K2 := H2fi — H2fi is of finite rank 0{h~^), has a compact support in the sense that K2 = 
X2-f^2X2 if X2 G Co° is equal to 1 on .6(0, R) for some sufficiently large R, and 

{H2,e - z)-^ = 0(1) : ^ ' — ^ D{H), uniformly for z G H. 

We repeat tlie proof of the Theorem 2 replacing by K2 and Kq{z) by K2{z) = K2{z — 

H2,e)~^ ■ Consequently S^log det ^1 + ^0(2:)^ is replaced by S^log det ^1 + K2{z)^ which is a 

non-holomorphic function. We treat this term as the term c^zlog det ^1 + Ki{z)^ in the proof 
of Theorem 2. □ 

Remark 8. The equation (6) shows that the spectral shift function ^{X,h) satisfies 

Im Wy^O 

(45) + [Ao,A]; /xGad(/i'i)}. 

In particular, for X G I\ad{Hi) the distribution ^(A, h) is continuous, and the function 

rj{X, h) - r?(Ao, h) = 4(A, h) - 4(Ao, h) - G [Aq, A]; /x G ad{Hi)} 

is real analytic in I. □ 

Repeating the argument used in the proof of [8, Theorem 4], the following theorem is a 
direct consequence of Theorem 2. 

Theorem 6. (Local trace formula) Let Q be an open complex, simply connected and rela- 
tively compact set satisfying the assumption (A^) such that / = J7 n M is an interval. 
We suppose that f is a holomorphic function in 17 and ifj G Cq°(R) satisfies 

' \ 1, d(/,A) <e, 
where e > and sufficiently small. Then 

tr[iijf)iH)]l= Yl fiz) + Ea,fAh), njith 
2eRes(i?i)nn 

\En,f,Ah)\ < M{il;,n)snp{\f{z)\; < didn,z) < 2e, lm{z) < 0}h~\ 
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7. Weyl Asymptotics 

In this section we obtain a Weyl-type asymptotic for the spectral shift function ^(A, h) 
associated to operators Hq and Hi = Hq + V. Here we assume that V is an electro-magnetic 
potential (8), 



ev. 



Hi = — aj{ichdj + eAj) + mc?P + i 

i=i 

In the following, we fix Iq C M\{±mc^} and choose Wq an open simply connected, relatively 
compact subset of $7 satisfying the assumption (A^) such that Iq = Wq fl R. 

For the /?,-pseudo-di£Ferential and functional calculus for the Dirac operator, we refer to 
([12], [10], [34], [19]). We recall that H^ = Op'f^{V^) and ip{H^,) are /?,-pseudo-di£Ferential opera- 
tors for a smooth function ip. The semi-classical symbol V^, is defined in (10). 

Let us introduce the intervals /i, I2 C Iq neighborhoods of Ai, A respectively such that, 
each A G /i U /2 is a noncritical energy level for H (see Definition 2). Let cpj G Cq°(M,M"'") be 
such that 

(46) (^1 = 1 on Ii, (/?2 = 1 on I2 and (fi+(f2 + f3 = '^on Iq. 

Consider a function e{t) G C^(] - Si,Si[), 9{0) = 1, e{-t) = e{t), so that the Fourier 
transform 6 oi 9 satisfies ^(A) > on M, and assume that there exist < eo < 1, Sq > 0, such 
that ^(A) > Sq > for |A| < cq. Next, we introduce 

(J^;7^^)(A) = {2Trh)-^ j (^*^^~^ e{t)dt = {2'Kh)-^e{-h-^X). 

To prove Theorem 3, we need the proposition: 
Proposition 10. For the trace involving H^, 



(47) 



tr 



0, 1, we have for A G /j. 



{J'-'9){X-H)<pj{H)\]=wj{X)h-^ + 0{h-'), j = l,2, 



with Wj{X) G C^{Ij) and 0{h ^) uniform with respect to X ^ Ij. 

Proof. Proposition 10 is closed to the calculation of the trace in [5, Section 4] and to the 
appendix of [8] for the Schrodinger operator. But, here we use a trick of Robert [10]. We fix 
J = 2 (it is similar for j = 1). The proof of (47) is obtained following these two steps: 

• First, we recall that X G I2 and Supp^(i) C [— 5i,(^i]. Let us write 



r = tr 



{J-^'e){X-H)MH] 



=tr 



27r/i 



,it(X-H.)h- 



<P2iH.)dt 



2^ / e'*^''"^(*)tr[e-'*^'^"V2(i?.)]o^i 



In the order to calculate the trace 

trifiHi) f{Ho)), for all / G C^iR\{±mc'}) 

we use [10, Proposition 3.2]. If we note W{h) = Q - li[Q,A{h)] with Q = Hf - H^, 
A{h) = l{x-hd:c + hd^ ■ x) and [Q,A{h)] = QA{h) - A{h)Q, we have 

(48) trifiHi) - fiHo)) = tT{W{h){Hf - m'^cY' fiHi)). 
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Applying the formula (48) for /(A) = e~'^*^^ V2(A), we have 

T = ^ y e'^^''~'e{t)ti (w{h){H^ - n?c^y^e-'^^'^~' ip2{Hi)) dt. 

Remark 9. Of course {Hl — m?c^)~^ is not well defined, however for f G (7o°(M\{±mc^}), we 
can define {H^ —m'^c^)~^ f{Hi) as the self-adjoint operator <f {Hi) where € Cq°(M) satisfies: 

(A2 - m2c^)-i/(A) for A / imc^, 
for A = ibmc^. 



^(A) 



• Now, we treat T following the analysis of [5, Section 4.2]. By /i-pseudo-differential cal- 
culus, we obtain the existence of a /i-pseudo-differential operator S which is trace class with 
symbol 

(49) s{x,y,^,h)eS\{xr\0~''). ViVGN, ,5>3, 

having compact support in ^ and in {x — y) (i.e. supp(j._j^)(s) = {x — y, 3 ^ ; {x,y,^,h) G 
supp(s)} is compact) and support in {(x,^); \x\ > R, (x,^) G T>^^(l2)}, with T>i the semi- 
classical symbol of Hi, so that 

r = ^tr (^J (^^^^~' e{t)e-'*"^^~' Sdt^ +0{h'^). 

Using Theorem 9 in Appendix A and the hypothesis of S by composition of Fourier integral 
operators, we obtain a Fourier integral operator lit = U^' +Uf , such that for \t\ < Si and Si 
sufficiently small, we have 

(50) \\Ut-e-''''^''-'S\\tr = 0{h'^), 

where the kernel of the operator / g**^'* ^9{t)Utdt is equal to K'^{x,y;h) + K~{x,y;h) with 

The amplitudes E^, satisfy 

E^{t,x,y,^;h) G S\{x)-' {O-""), VAT G N 
and is compactly supported in ^ and in {x — y). 

Using the Taylor sum formula for the functions ^>^(t, x,^) in a neighborhood of i = 0, we 
have: 

$±(t, X, = a; • ^ - tHfix, + O(t'). 
We will deduce that T = T+ +T-, with 

(2^ /// e'^*^'^**^*''''^^~'"^^''"'^(*)^^(*'^'^''^;^)^*^^^^ + o(^°°)- 

Moreover, the symbol E^{t,x,x,^;h) has support in {(x,^); |x| > R, |^| < Ci, {x,^) G 
T^i^ih)}, so that for all a and \t\ < Si, we have 

(51) \d''E^{t.x.x.^;h)\ <Ca{x)-^, S>3. 

The last estimate enables us to calculate T by using an infinite partition of unity 

^ ^-(0; - a) = 1, Va; G M^, 
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where * G Cq°{K), ^ >0, K being a neighborhood of the unit cube. Consequently, for every 
fixed h E]0,ho], we have 



{2Trhf rn^^ J J J 



X V ^{x -a)E^{t,x,x,C;h)dtdxdC + 0{h°°) = Urn l^ + 0{h°°), 

|a|<m 

and we reduce the problem to the analysis of the integrals I^. Concerning the phase function, 
we observe that 

(52) tX + X, e) - X • e = t(A - iff (x, + 0{t)), 

where 0{t) and its derivatives are uniformly bounded on the support of d{t)E^{t,x,x,^; h) 
since the derivatives of {^^{t,x,^) — x ■ ^) are bounded on this set. 

Now we look for critical points of the phase function {t\ + x, ^) — x • ^). Putting the 

derivative with respect to t equal to 0, we see that Hf{x, ^) = \+0{t). Since dx^fHf{x,S^) ^ 0, 
when H^{x,^) = A, and putting the derivative of the phase function t{X — + 0{t)) 

with respect to H^(x,(,) equal to 0, we have 

t = 0(^2). 

Then the phase is critical for \t\ small precisely when t = 0, A = H^. Near any such critical 
point we choose local coordinates t, H^{x,^), wi,--- and consider the Hessian of (52) 
with respect to t, Hf{x,^) at the critical point: 

★ -1 
-1 

This is a non-degenerate matrix of determinant —1 and of signature 0. By the stationary 
phase method we obtain 

lt. = j^ I , E ^{x-a)E^{Q,x,^,\-h)L^{dw)+0{h% 

where L^{dw) is the Liouville measure on A = Hf^ and the remainder 0{h~'^) is uniform with 
respect to A G /2 and m G N. Here V'^(A) G C'q°(/2)- Taking the limit limm-^oo obtain 

an asymptotics of T. □ 

Lemma 9. With the above definitions of 9{t), ^(A, /i), <^j(A), Ij, j = 1,2, we have 

(53) / J'f^'e * ifj^'i,,, h)dii - t ^j{iJi)i'{n)dii = 0{h-^), A G 

Proof. We deal only with the analysis of (53) for j = 2 since that of j = 1 is similar. 
According to Theorem 2, there exists a holomorphic function r{z,h) in O such that for all 
A G /o = W^o n M, we have 

ei\,h) = -lmriX,h)+ Yl ,'^""'',2 + E -^-W' 

TT ' 7r|A — ^-^ 

weRes{Hi)nn «;eRes(J/i)n/o 

Imw^o 
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where r{z, h) satisfies the following estimate: 

(54) \r{z,h)\<C{W)h-^, zeW, 

with CiW) > independent of h. Let us denote 



G^2W = - Im r{n,h)ip2{n)dii, 

TT J-oo 



weRes{Hi)nn '°° «;eRes(ifi)n]co,A] 
Imuj^o 

Using the Cauchy inequalities and (54), it follows easily that 

G;^(A) = Oih-') and G';^(A) = ©(/i'^), 
and we immediately obtain 

(56) J^t^'e*G'^^-G'^^ = 0{h-'). 

Now, we want to apply a Tauberian theorem (see [34, Theorem V-13 ]) for the increasing 
function Mtp2(X). For this purpose, we need the estimates 

(57) M^,(A) = 0(/i-3), ^(^-i^*M^J(A) = 0(/i-3), VA G R, 
and the equality M^^{fj,) = G^^{iJ,) = 0, fj, < infl2- 

The first estimate in (57) follows easily from the equation (55) with the upper bound of the 
number of the resonances in O, (see Theorem 1), and the second follows from (47) and the 
equation 

^{^h'(^ * M^J(A) = * ^2e'(A) - ^i^h'() * G^J(A)- 
Then, according to the Tauberian theorem we have 

{J'-^e * M^,)(A) = M^,(A) + 0{h-^), 

this enabled us to obtain 

r ip2m\fi)dfi = M^,(A)+ r G'^^ifi)dfi 

J —oo J —oo 

= t ^iJ'h'9*M^,+J^-'e*G^,)in)dfi + 0{h-^) 

J —oo 

= r :F^^e*ip2eii^,h)dii + oih-^). 

J —oo 

□ 
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Proof of Theorem 3. For Ai G /i, A G 72, using the functions defined in (46), we have 

/A pXi 
-oo J —oo 

f Ai t-X 

^i{lAi'{l^^h)dii+ 1 ip2{n)^' in, h)dfi 

-oo J —oo 

r-A 

+ 



/ 'P3{f^)C'{l^,h)di2. 



Since (pj = on I^-j for j = 1,2, the first term (resp. the second term) is indepen- 
dent of A G /2 (resp. Ai G h) and is equal to tr[(pi{H,)]Q = C{ipi)h~^ + 0{h~'^) (resp. 
iv[<f2{H)W = C{(p2)h~^ + 0(/i~^)), where C{ipj) is a constant depending on ipj for j = 1,2. 
Since = on Ij, j = 1, 2, the last term is independent of A G I2, Ai G h and is equal to 
C{(ps)h~^ + 0{h~'^), where C{(ps) is a constant depending on ip^- The proof of this results is 
based on the functional calculus in the framework of /i-pseudo-differential operators. 

Using the equations (47), (53) and (58) we complete the proof of the asymptotic expansion 
(13) by writing 



{^H'0^{^,m>^) = ((^-ie)(A-.)^,(-),0 



= tr 





-3 , Al/j,-2 



(59) = WjiX)h--' + Oih-^), j = 1,2. 
It remains to compute the Weyl term (14). 

According to the definition of the spectral shift function C(A, /i) in (5), we have: 

(60) (C'(A, h),<p{X)) = tv{^{Hi) - ip{Ho)), ^{X) G 



We use weak asymptotics which is a direct consequence of functional calculus in the framework 
of /i-pseudo-differential operators, as settled in [12], [34], [10]. We find 

= Op%{V^), = 0, 1, (Pz. defined in (10)), 

and 

tv{<f{H,)-^{Ho)) = h-'J2^j{ip)h^ 

= hSo{^) + 0{h-^), 

with 70(99) = {2tt)-'' /jj,,/g3 tr (ip{Vi{x,C)) - (p{Vo{x,C))) dxd^ 
(tr(^) is the trace of the matrix A). 

The matrix 'Di^{x,^) is Hermitian and has two eigenvalues H^{x,^) (see (11)), then 

tr (^(Pi) - (^(Po)) = 2 {ip{Ht) + ^{H^) - v{H+) - ^{H^)) . 
According to the asymptotic expansions (13) and (60) we obtain 

u;(A, Ai) = ■u;(A) - ■u;(Ai), 
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with 



MA) = -^/(/ dC- [ dC- [ di+ j dAdx. 

471" J'8?\Jh+(x,(,)<\ JH+(x,i)<\ •JH:[{x,0>\ JHq{x,^)>\ J 
Putting = — veA{x) for = 0, 1 and ^ = ra; (a; G 5^), we get 

<±\ <^ Kiz + l^ic +v ) 1 ±v < ±A 



1 



thus 



Jh+{x,o<\ Jh-{x,o>\ V ^ ^ / + 



3 \ 2 2 



+ 



and 



-/ I d^ = T^(A2-(mc2)2)i for ±A>0, 

with = max(x, 0) and = max(— x, 0) for x G R. □ 

Remark 10. Theorem 3 can he extended to the pair operators {Hi = Hq + Vi, H2 = HQ + V2), 
where the potentials Fi, V2 are electro-magnetic potentials 

satisfying the assumption (Av) with 6 > (or \\Vj{x)\\ — > 0) and the potential V = V2 — Vi 
satisfies the assumption (Av) with 6 > 3: 

For all A, Ai noncritical energy levels for Hi, H2 such that ±mc^ ]Ai, A[ and h €]0, /io[, we 
have the asymptotic expansion 

(61) e(A, h) - e(Ai, h) = w{\ \i)h-^ + o{h-^). 

Here the 0{h^^) is uniform for A (resp.Xi) in a small interval I2 (resp. Ii). The first term 
u;(A, Ai) G C°°(/2 X Ii) is given by 

w{X, Ai) = w{X) — w{Xi) 

with, 

(62) '^^^^ = ^ [_[W+{X,V^,V_)-W+iX,V^,V_)]ldx 



3 

where W± (A, a, h) = ((x - ^^^) ^ - (mc^ + 



In this setting, we does not have a formula like (48)- But it could be possible to use the 
approach to Bruneau-Petkov in [8]. For that we need more informations on the approximation 
of the propagator e~^^^^^ by the Fourier integral operator Uf 
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8. Resonances in small domains 

In this section, we assume that the hamiltonian H = Hq+V, where V is an electro-magnetic 
potential: 

3 

H = — aj{ichdj + eAj) + mc^P + ev. 

8.1. Upper bound for the number of resonances in domains of width h. We adapt, 
for the Dirac operator. Theorem 1 of [9] which rests on a representation formula for the spectral 
shift function (see Theorem 2). 

Theorem 7. Suppose that each A G [Eo,Ei] is a non-critical energy level for H. Then for 

h g]0,/io], the following assertions are equivalents: 

(i) There exist positive constants B, C, b, ho, such that for any A G [Eq — b,E + b], 
h g]0, ^o] o.nd h/B < p < B , we have 

#{z eC: ze Res{H), \z - X\ < p} < Cph~^. 

(a) There exist positive constants Bi,Ci,ei,hi, such that for any A G [Eq — Ei, Ei + ei], 
h g]0, ^i] and h/Bi < p < Bi, we have 

\^{x+p,h)-a^-p,h)\<cph-^. 

As a consequence of Theorem 7, we have an upper bound 0{h~'^) for the number of the 
resonances for the semi-classical Dirac operator close to a non-critical energy level in a domain 
of width h: 

Proposition 11. Assiime that V is the electro-magnetic potential (8) satisfying the assump- 
tion (Av) with 5 > 3. We suppose also that each A G [Eo,Ei] is a non-critical energy level for 
H. There are positive constants C, B, b, ho such that for any A G [Eq — b, E -\- b], h g]0, ho] 
and h/B < p < B, we have 

#{zeC: ze Res{H), \z - X\ < p} < Cph~^. 

Proof. It follows from Theorem 7 and equation (13). □ 

8.2. Breit-Wigner approximation. In this part, we consider small domains of width h, and 
we prove a Breit-Wigner approximation for ^{X,h) (see [31], [32], [16], [6], [8]). Let ri{X,h) be 
the real analytic function defined by 

r?(A, h) = e(A, h) - G [Eo, A] : /X G ad{H)}. 

Using Proposition 11 and the arguments used in [8, Section 6], we obtain a Breit-Wigner 
approximation for the derivative of the spectral shift function ^(A, h). 

Theorem 8. (Breit-Wigner) Assume that V is an electro-magnetic potential (8), for any 
X G [Eo,Ei] non-critical energy level for H, < p < h/B, < Bi < B, and h sufficiently 
small, we have 

r]{X + p,h)-r]{X-p,h)= Yl '^C-{yJ,[X- p,X + p]) + 0{p)h-^, 

weRes{H) 
Imw^O, \w-\\<h/Bi 
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where B > is a positive constant and loc_ is the harmonic measure 
Loc_(w,E) = -- [ ^'^^'^i dt, EcR = dC-. 

T^JE\t-w\ 

Using Theorem 7 and repeating with httle modifications the arguments used in [7, Section 
6], we obtain the following corollary which implies also trace formula in small domains. 

Corollary 1. Under the assumptions of Theorem 8 and supposing that [Eq^Ei] is a non- 
critical energy level for H, for each E G [Eq,Ei\ there exist constants C2 > Ci > 0, /iq > so 
that for |A — £"1 < Ci/i, h g]0, /iq], we have 

(63) a\h) = -\ E u>^)+o{h-'). 

weRes{H) wGcTaiH) 
\E-w\<C2h \E-w\<Cih 

Here 6w{-) is the Dirac mass at w gM.. 

Appendix A. Construction of Ut 

In this appendix, we construct a parametrix at small time of the propagator of the Dirac 
equation in an external electro-magnetic field 

ihdtip = Hit/j, 

with Hi = Hq + V . Here Hq is the selfadjoint operator defined in (1) and V is an electro- 
magnetic potential (8). 

Theorem 9. (Approximation of the propagator) There exist 5\ > small enough and 
a Fourier integral operator lit = + Ut '^^^^ 

^^^^^^ = (2^ / j 

defined for \t\ < 5i such that: 

• The amplitude E^{t,x,y,C;h) £ cS°(l). 

• \\Ut - e-^*"^^'''\\ = 0{hP^), uniformly for \t\ < 61. 

• The phase function $^(t, x,^) — x ■ ^ and their derivatives dfdxdj x,$,) — x ■ ^) 
are uniformly bounded for {t,x,$,) € [— 5i,(5i] X X B{0,Ci), (a,/?, 7) 7^ (0,0,0) and 
Ci > (see(68)). 

With a different approach, such result has been obtained by Yajima [46] for a scalar electric 
potential (f+ = V-). 



(64) 



Proof. We consider the equivalent problem for Ut 

ihdtUt-HiUt =0, 



We solve this problem using the B.K.W. method. We assume that the kernel of the operator 
Ut is Kt, where 
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with E{t, X, y, ^; h) = Eo{t, x, y, C) + hEi {t, x,y,0 + ---- 

Thus, if we look for E{t, x,y,$,; h) having the asymptotic expansion above, it is enough 
solve (in some fixed neighborhood of t = 0) the sequence of equations 



(65) { 



= [dt^{t, X, + ca • - ea- A + mc^(3 + ev) Eq, 

i{dt + ca ■ Vx)Ej =(5t$(t, x, ^) + ca ■ Va;$ - ea- A + mc^P + ev) Ej+i, 

Eo{0,x,0 =h, 

Ej{0,x,0 =0, forj>l. 



(66) 



(67) 



On the support of Eq, we deduce the eikonal equation 

det {dt^t, X, + ca ■ -ea- A + + ev) = 0, 
^O,x,0=x-^. 

The system (66) is equivalent to 

dt^^it, x, + Hfix, V.$) = 0, (see (11)), 

<i>±(o,x,e) = .x-e. 

The last system can be solved by Hamilton- Jacobi method (see [2]) and all derivatives 

(68) d?d^dJ{^Ht,x,0-x-^) 

are uniformly bounded for (t, x, ^) G [—(^1,(^1] xR^ X 5(0, Ci) and (a,/3,7) / (0,0,0). 
Using the Taylor formula in a neighborhood of t = 0, the two solutions of (67) satisfy: 

$±(i, x,0=x-^- tHtix, + o{e). 
Then lit = + , and the kernel of the operator lit is Kt = + , with 

We look for the amplitude E^{t,x,y,^;h) having an asymptotic expansions according to t 
power of h: 

E^{t, X, y, + hEf{t, x,y,0 + -- - ■ 
Consequently, the coefficients E^{t,x,y,^) are the solutions of the transport equations 

= (af<I>^ + ca ■ Vx^^ - ea- A + mc^fS + ev) E^, 
i{dt + ca - Vx)Ef = + ca - Vc,$^ -ea-A + mc^P + ev) Ef^^, 

E+{0,x,0 + E-{0,x,0 =0 fori>l, 

E^{o,x,o =nf(x,o, 

with Uf{x,0 defined by (12). 



(69) { 



Resolution of (69). 

Let us denote hy L = dt + ca - Vx, with a - Vx = Yl^=i (^jdxj • The matrix 

= dt^^ + ca ■ Vx^^ - ea - A + m.c^ (3 + ev, 
is Hermitian and has two real eigenvalues linearly independent with multiplicity 2. 
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First, we multiply the system (69) by the column- vector A^i = (1,0,0,0)^, the subscript f 
designates the complex conjugate of the transpose. We denote 

(70) Ef^^ = EfNi for j = 1, 2, • • • , ^^^^^(O, x, = nf (x, 

Since det(A^^) = 0, there exist and r^, left and right eigenvectors of the matrix M^, 
corresponding to the eigenvalue zero, such that 

(71) M±r± = 0, l^M^ = 0, Z^ = (r±)t, A; = 1,2, 
(here is a column- vector and is a row- vector). We choose 
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(72) Z±r± = {T2pfu^)6,k, u,k = l, 2. 

Here , and are defined by 



pi, 



w^ = ±pf + ip^, wt 



ipf, 



where = m(? + 



Pb 



of = cd^,^^ - eA^, for 



1,2,3. 



It is easy to see that the vector-valued functions r^(t, x,^) and l^{t,x,^) can be chosen to 
be smooth in t and x and different from everywhere. All the derivatives of r^,l^, k = 1,2, 
are uniformly bounded for {t,x,$,) G [— ^i,<5i] X X 5(0, Ci). Then it follows from the first 
equation in (69) that 

where 0-^1,0-^2, ^^^'^ scalar- valued functions. If we multiply the second equation in (69) for 
_7 = on the left by for A; = 1, 2, we deduce the following differential equations for a^^: 



ltL{<irT) + ltL{a^or^) = 0. 



We conclude 



(73) 



+ZfL(r±)4i + /±L(rf)ao±2 = 0, 

ifrfdtiai^) + cE5=i l2(^jridxj{(^o,2) + cE5=i ^f^j^f ^x,((tJi) 

+ltL{rt)ai, + liL{rf)ai, = 0. 



We now use Lemma 10 (see below) in the system (73). Since pf ^ 0, = pf ^ pf ^ then, 
after multiplying (73) by {'^2pfu^)~^ , the system (73) can be written as 



(74) 



D^ai = M^ai := (T2p±u±)-i 



^0,1 
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with = dt + a^ ■V^ = dt + EjLi afit, x)d^., and 

Thus the function a"^^ can be found if its value is known for t = 0, and it is as smooth as 
C7^j^.(0, X, ^) (for more details, see a method for solving a similar equation in [37]). The equality 



<i (0, X, = (^ti X, e)rf (0, X, + (^±2 (0, X, 0^2^ (0, X, = nf (x, OiVi , 

gives the value of cr^ at t = 0. 



Since, all the derivatives of ctq ^, , for k = 1, 2, are uniformly bounded, then all the deriva- 
tives (dfd^djE^^^) are uniformly bounded for (a, ^,7) e N x x N^. 

It follows from the second equation in (69) for j = 0, that 

i.e., E^^i = (J^^ir^ + (^^^^2^2 + ' where af^ is a scalar-valued function for A; = 1,2, and 
is expressed in term of LE^^. To find cr^^ it is sufficient to multiply the second equation in 
(69) for j = 1 on the left by for = 1, 2. Then 



< 



[ lfL{at,rt) + lfL{at2rf) + ltL{hf) = 0. 



From this equation, o"^^ can be found if the function o"i yt(0, x) is known. With the same 
process, for all ,7 = 1, 2, • • • , we obtain 

4i'-t + 42^t + hf =Ef^„ 
lfL{af^^rf) + ltL{ayt)+lfL(hf) =0, 
[ ltL{af^^rf) + ltL{cjyt)+ltL{hf) =0. 

For t = 0, _7 = 1, 2, ■ • • , we have 

and the quantity is determined if E^^, Ef^, ■ ■ ■ , E^_^ are known. Solving the differential 
equation for (^f = ^ , we find these functions for all sufficiently small t. 

Repeating this block of calculus, multiplying by N2 = (0,1,0,0)^, A/3 = (0,0,1,0)^ and 
Ni = (0,0,0,l)t instead of A^i in (70), we find Ef^^= EfN2, Ef^^= EfN^ and Ef^^= EfN^. 
Consequently, we have: 

Proposition 12. There exists a family of matrices 

= (^ti^ ^^2, Ef^s, Ef^d, forj>0, 

solution of (69). Moreover, for all j > 0, E^ G C°° and all derivatives {d^dxdjE^) are 
uniformly hounded for all (t,x,^) G [-61, 5i] X M3 X B{0,Ci) and (a,/3,7) G N x x N^. 
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Consequently, Borel process provides a symbol E^{t,x,y,^;h) G S'^{1) with compact sup- 
port on ^ and {x — y) with E^{t, x, y, ^) + hEfit, x,y,^) + ■ ■ ■ its asymptotic expansion. 

Desired estimate. 

Next, we remark that for all TV G N: 

X T.U^^hL{Ef)+M^Ef)h^ 
(75) = PN{t,x,i-h)h'' , 

and all derivatives D'^ ^Pj^{t,x^^; h) are bounded as /i ^ for all a. Then for all A'' G N, 

{ Uo =I + 0{h^), 

thus 



^^^^ \ Wo =I + 0{h% 

where 0{h^) is uniform on t and corresponds to the norm in £(^2). Then we get: 

(78) \\Ut-e-'^"'^~'\\=0{h'^). 

□ 

Lemma 10. Under the notations used above, we have 

(79) lfajrf = l^ajr^, lfajrf = lfajr^ = 0, j = 1,2,3. 

Proof. As Rubinow and Keller in [37] let us work in a general situation. 
We consider the n Hermitian matrices and n real scalars p^, fj, = 1, • • • ,n. Let G be the 
Hermitian matrix defined by 

n 

Let A be a multiple eigenvalue of G and -Bi, • • • ,Bq, a set of corresponding orthormal eigen- 
vectors which are differentiable functions of . Then 

(80) B}Bk = 5jk, 

(81) GBk = XBk. 

The subscript f designates the complex conjugate of the transpose. 

If A(pi, • • • ,Pfj,) is differentiable, we differentiate (81) with respect to and obtain 

(82) M,B, + G^ = -^B, + xf^. 

The multiplication of (82) on the left by Bj, the use of (81), and the fact that G is Hermitian 
yield 

(83) B}M,Bk = g^^Sjk. 
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In order to treat our case, we take 

5 

where Mj = aj for j = 1,2,3, M4 = /? and M5 = I/i are Hermitian matrices {o(j,P are the 
Dirac matrices) and are five real scalars. 

We also take = pf ± ^ (pf)'^ + (pf)^ + (pf)^ + (Pa)^ -^^ the point which have the 
coordinate pf = cd^.^^-eAj for j = 1,2,3, pf = mc^ + "^^+~^-\ pf = + . 

When satisfies (66) and (67), r^,r^ are two orthogonal eigenvectors of corre- 
sponding to the eigenvalue = X^{F^) = 0. Since |e(f+ — V-)\ < 2mc^ (see (9)), is 
differentiable near the point F^. Now, we apply (83) with B^- = if and = . After the 
normalization of r. , / • we obtain 



dp^ 



^^{^2p^u^)6jk, 



and we get the lemma. □ 
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